Macroscopic quantum effects generated by the acoustic wave in a molecular magnet 
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We have shown that the size of the magnetization step due to resonant spin tunneling in a 
molecular magnet can be strongly affected by sound. The transverse acoustic wave can also generate 
macroscopic quantum beats of the magnetization during the field sweep. 
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Single-molecule magnets(SMMs) have attracted much 
interest because they provide possibility to observe quan- 
tum effects at the macroscopic scale. Among these effects 
are step-wise magnetization curve caused by resonant 
spin tunneling [lL |2|, topological interference of tunnel- 
ing trajectories [3|, and crossover between classical and 
quantum superparamagnetism 0, Q ■ Landau-Zener the- 
ory has been used to describe spin transitions that occur 
during the field sweep @. It has been recognized that 
spin-phonon interactions play an important role in the 
dynamics of spins in molecular magnets 0, @] ■ Possibil- 
ity of Rabi oscillations of spins caused by the acoustic 
wave has been studied In recent years the effect of 
sound on molecular magnets has been explored in ex- 
periment (lo| . In this Letter we show that sound can 
significantly affect the size of the magnetization step due 
to resonant spin tunneling. In the presence of the field 
sweep an acoustic wave can also generate quantum beats 
of the magnetization of a macroscopic sample. We com- 
pute the parameters of the sound that are necessary to 
observe these effects. 

The effect we are after is illustrated in Fig. 1. The 
lower curve shows how the sound modulates the distance 
between spin energy levels in the absence of the field 
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FIG. 1: Schematic diagram of the time dependence of the 
distance between spin energy levels. Thin solid line: The 
effect of the acoustic wave without the field sweep. Thick 
solid line: Field sweep is modulated by the acoustic wave. 



sweep. In this case the phase of Rabi oscillations caused 
by a propagating sound wave depends on coordinates 
such that the oscillations average out over the volume 
of the sample if the latter is large compared to the wave- 
length of sound. In the presence of the field sweep (the 
upper curve) the phase of the Rabi oscillations is still 
a function of coordinates. However, the Landau-Zener 
probability of spin transitions that contribute to the os- 
cillations depends on the rate of the field sweep. That 
rate becomes modulated by the sound. Consequently 
the regions of the sample that contribute most to the 
dynamics of the magnetization add their contributions 
constructively. The resulting oscillations of the magnetic 
moment of the sample can be observed in a macroscopic 
experiment. 

We consider a crystal of single-molecule magnets with 
the Hamiltonian, 



^smm = —DS Z — gnBH z S z + Tit 



(1) 



where Si are Cartesian components of the spin operator 
and D is the second-order anisotropy constant. The sec- 
ond term is the Zeeman energy due to the longitudinal 
field H Zl with g being the gyromagnetic factor and [Ib 
being the Bohr magneton. The last term includes the 
transverse magnetic field and the transverse anisotropy, 
which produce level splitting. Local rotation produced 
by a transverse acoustic wave of frequency uj = Ctk, wave 
vector k, and amplitude Uo, polarized along the y axis 
and running along the x axis, is given by 11 1 



S(j)(r) = —kuo cos(kx — u>t)z . 



(2) 



Due to the rotation of the local anisotropy axis by sound, 
the spin Hamiltonian becomes [8| 



H = e 



rtSMMe 



(3) 



The simplest solution of the problem for an individual 
spin can be obtained in the coordinate frame that is 
rigidly coupled to the local crystallographic axes. The 
wave functions in the laboratory and lattice frames, 
and |*( lat )), are related through 



(lat)\ = jStfr-S 



(4) 



2 



A 

N 

CO 

V 



6.4 
6.0 
5.6 
5.2 
4.8 
4.4 



'\^. i i^ l 'ii l iCl , '...,''; , .!'KI |,i, it^w« J i;! l 



-p=0.9,q=0.247 
p=1.0,q=0.2 
-p=1.1,q=0.165 



30 



35 



40 



45 
X 



50 



55 



60 



FIG. 2: Color online: Time dependence of the average spin 
of the sample for given values of p and q at S — 10, M = 0, 
and 7 = 0.02. Inset: Magnetization step and oscillations in 
the wider range of r. 
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FIG. 3: Color online: Final magnetization vs p for the M = 
step with and without acoustic wave at different 7 and q. 
Inset: <S Z ) vs r at 7 = 0.01, <? = 0.2 and p = 0.9. 



while the spin Hamiltonian in the lattice- frame is given 

by HHE3 



w (lat) = ^ 



SMM 



as • n , 



with 



f2 = ^(/> = - — uq sin(kx — ui)z . 
2c t 



(5) 



(6) 



We are going to solve the problem locally for each spin 
in the lattice frame and then use the above formulas to 
obtain the solution for the entire crystal in the laboratory 
frame. 

In the absence of transverse terms the energy levels of 
the Hamiltonian ([TJ) are 



-Dm —g^BH z m. 



(7) 



where S z \m) = m\m). Close to the resonance between | — 
S) and \S—M) the Hamiltonian ((5]) can be projected onto 
these states, resulting effectively in a two-level model: 



H. 



(lat) 



1 



s =--A& x -5E(a z +I), 



(8) 



where 



g\iBct -I ^— sm(te — wt) 



a z = \S-M)(S-M\-\-S)(-S\ 
a x = \S-M)(-S\ + \-S)(S-M\ 
f=\S- M)(S - M\ + I - 5) (-5| , 



(9) 



A is the splitting of the resonant levels, c = dH z /dt is 
the field sweep rate, and 



2c t 



(10) 



is the Rabi frequency. Treating x as a parameter, we 
express the corresponding wave function as 

\^\t)) = b S - M {t)\S M) + b_ s (t)\ - S). (11) 

and solve the time-dependent Schrodinger equation, 
r (lat). 



.(lat), T (lat) 



at ~ off 1 off 



(*)), 



(12) 



that at M = becomes equivalent to the following two 
coupled differential equations: 



dr 
dbs 



2iS 
'bs, 



qp 



jt — sin(pr — kx) 

o 



dr 2' 
where we introduced dimcnsionlcss 



(13) 



r = t 



7 : 



A 2 ' 



P = 



A"' 



q =^. (14) 



We consider samples of length that is large compared to 
the wavelength of the sound. The expectation value of 
the z-projection of the spin at M = is given by 



^(t)\S z \^(t)) = S\bs(r)\ 2 + (-S)\b^ s (r)\ 2 



(15) 



Eqs. (fT3|) have been solved numerically. Fig. [2] illus- 
trates situation when the field was changing at a constant 
rate 7 and a pulse of sound was introduced shortly be- 
fore reaching the resonance between the | — S) and |5) 
states. The tunnel splitting is assumed to be sufficient 
to produce transitions between these two states. The 
most striking feature of the magnetization dynamics ob- 
served in simulations are the beats which are in line with 
the idea outlined in the introduction. Fig. [3] shows the 
p-dependence of the final magnetization on crossing the 
step for various values of 7 and q. This strong dependence 
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FIG. 4: Color online: (S z ) vs r for rj 2 = 0, 0.1, 0.5, and 1.0 
at S = 10, p = 0.9, q = 0.247, 7 = 0.02, and M = 2. Inset: 
Oscillations between r = 20 and r = 50 at 772 =0 and 0.1. 
Note that the beats and the oscillatory behavior disappear as 
772 increases. 



of the magnetization step on frequency and amplitude of 
the acoustic wave, as well as on the sweep rate, is one 
of our main results. We believe that it should not be 
difficult to observe this effect in experiment. 

Another possible experimental situation corresponds 
to the sample initially saturated in the | — S) state, after 
which the acoustic power of frequency u> ~ A/h is applied 
to the crystal and maintained during the sweep. Here 
\S — M) is the level that at a given sweep rate provides 
significant probability of the transition when it is crossed 
by the | — S) level. In order to study such a problem, 
we need to know the rate of relaxation of the \S — M) 
state to the lower energy states. Defining r s-m+i,S—m 
as the rate of the \S — M) — > \S — M + 1) transition 
and introducing t/m = K^s-m+i,s-m/ A, we obtain two 
coupled differential equations: 

dbs-M . , i, 

= izbs-M + 2 b -s, 



(It 
dr 



(16) 



where 
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FIG. 5: Color online: (S z ) vs r for M = 0, 1, and 2, at r] M = 0, 
S = 10, p = 0.9, q = 0.247, and 7 = 0.02. 



in Fig. |4] As the damping increases, the magnetiza- 
tion jump becomes more pronounced but the oscillatory 
dynamics disappears. The comparative behavior of dif- 
ferent resonances is shown in Fig. O 

Let us now study the optimal conditions for the ob- 
servation of the macroscopic acoustic Rabi effect studied 
above. Defining u) q (r, x) = S[jT—(qp/S)sm(pT—2Trx/A)] 
where A = 2n/k is the wavelength of the sound, the cou- 
pled differential equations (I13|) can be written as 



2iui„bs — [ 2i(j„ — — ] bs = 



(19) 



where b$ = dbs/dr, bs — d 2 bs/dr 2 and so on. Introduc- 
ing 6s = ds exp[w(r)] and selecting v = u> q , we get 



d s + [ -iuJ q + uj q + -) d s = ■ 



(20) 



which describes damped oscillations. At q ^ we have, 
e.g., for x = and x — A/2 



w 9 (t,0) = S 



7T 



qp 



sin(pr) 



u) q (r, A/2) =S 1 T +1T sin(pr) 



qp 



(21) 
(22) 



respectively. This implies that each frequency gener- 
ates slow and fast oscillatory regions due to the sinu- 
soidal function, and they show different damped oscil- 
latory structures in a given range of r. In other words 
if £Jg(r, 0) is larger than lu = Sjt in some range of r, 
u) q (r, A/2) is smaller than u> , and vice versa (see Fig. 
[I}. At r ~ the frequencies are approximately given by 
(S z ) = -2S\b_ s (r)\ 2 -M\b S -u{T)\ 2 + S (18) u 9 (t,0) ~ S( 1 -qp 2 /S)T and w,(t, A/2) ~ S( 1 +qp 2 /S)T 

Introducing tan#i = u> q (r, 0), tan#2 = w g (r, A/2), and 
w = qp 2 1 'S, we get 



z=(2S- M) [jt - (qp/S) sin (pr - kx)] +ir] M /2. (17) 

As the lifetimes of the excited states with \S — M + 
1),...,|5 — 1) are shorter than the lifetime of \S ~ M), 
their contributions to the above equations can be ne- 
glected. Then 



We solve Eqs.(|T6|) numerically for selected values of 7 
and t)m- In the overdamped case, Ts-m+i,S-m ^> A, 
we find no Rabi oscillations. For the underdamped case, 
Ts—M+i,S—M ^ t ne numerical solution is illustrated 



|tan(0i-0 2 )| = 



2w 



I + 7 2 



(23) 
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FIG. 6: Time dependence of the magnetization of the sample 
for M = at S = 10, 7 = 0.006, q = 0.01, and p = 2.45. 



which increases monotonically in the range of < w < 
1 +7 2 . Under this condition, let us first consider two lim- 
iting cases: »»7 and w <g; 7. In the first case {0\ — 9 2 ) 
increases with w, and thereby \u) q {r, 0) — w q (r, A/2) | also 
increases, which is a less favorable situation for the beats. 
In the second case we have B\ ~ 62, which results in 
u)q(r, 0) ~ w g (r, A/2) ~ This also is not a favorable 
situation for the beats because it does not generate slow 
and fast oscillatory regions for x = and x = A/2, as dis- 
cussed previously. The optimal condition for pronounced 
beats is then 



7 



S 



(24) 



We shall now discuss what the above condition means 
for experiment. It is easy to see that Eq. (|24[) is equiva- 
lent to 



Uo 

A 



g 

ttS 



(25) 



The validity of the continuous elastic theory that we em- 
ployed requires uq <C A, that is, one needs to satisfy the 
condition q < 1. This is not sufficient, though. Since 
experiments on molecular magnets require temperature 
in the kelvin range or lower, one should also be con- 
cerned with the power of the sound. It should be suffi- 
ciently low to avoid the unwanted heating of the sample. 
The power per cross-sectional area of the sample is given 
by P/A = ^puQUJ 2 Ct- For, e.g., the parameters of Fe-8 
molecular magnet we find that the optimal conditions of 
the experiment require sound of frequency / = 0.5MHz 
-1MHz and power in the range 100W/cm 2 -200W/cm 2 in- 
troduced into the sample simultaneously with the field 
sweep of lkG/s. Time dependence of the magnetization 
under these conditions is shown in Fig. El 



In conclusion, we have demonstrated that the size of 
the magnetization step due to resonant spin tunneling 
in molecular magnets can be strongly affected by sound. 
The acoustic wave can also generate macroscopic quan- 
tum beats of the magnetization during a field sweep. The 
required frequency (MHz) and power (O.lkW/cm 2 ) of the 
sound, and the required sweep rate (lkG/s) are within 
experimental reach. 
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